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Abstract 



A reciprocal space formulation of the Cluster Variation is used in order to 
extract effective pair interactions in alloys from experimental short-range or- 
der diffuse intensities. The method is applied to the analysis of the short 
range order contribution to the neutron diffuse scattering of a Fe — 19.5%A/ 
single crystal. A detailed comparison with real space methods is carried out 
for three different levels of approximations. For the highest level of approxi- 
^ . mation used in this study, effective pair interactions up to fifth neighbors are 

H ! obtained. 
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I. INTRODUCTION 

The diffuse scattering of X-rays and neutrons remains tfie tecfinique of clioice for quan- 
titative studies of short-range order (SRO) in alloys. In particular, the part of the diffuse 
intensity related to substitutional SRO has played a prominent role in alloy theory. This 
SRO intensity contains important information on the effective interactions between atoms 
which, in turn, can be used to characterize thermodynamic properties of the alloy, such as 
relative phase stability, phase diagrams and antiphase boundary energies. 

Experimental measurements of SRO intensity have evolved significantly from the studies 
in polycrystalline materials of the early 1950's. At present, accurate measurements are car- 
ried out in single crystals at high and low temperature. For example, recent work has focused 
on the use of high flux neutron and tunable synchrotron X-ray sources.aLl On the theoretical 
side, the study of SRO intensity in alloys and its characterization through diffuse scattering 
dates back to the pioneering work of Cowley,@ KrivoglazB and Clapp and Moss.MliiJ Out of 
these studies emerged what is presently known as the Krivoglaz-Clapp-Moss (KCM) formula. 
This formula establishes a simple relationship between the experimentally observed diffuse 
scattering and the Fourier transform of the effective pair interactions. The KCM formula 
corresponds to a treatment of concentration fluctuations at the lowest level of approximation 
for the alloy configurational free energy, namely the single site mean-field approximation, 
also known as the Bragg- Williamsc3 and the molecular field approximation. 

More recently, effective pair interactions have been obtained from experimentally de- 
termined SRO intensities using higher approximations based on Monte Carlo simulations 
and the Cluster Variation Method (CVM).li2l The first calculation of SRO intensity carried 
out directly in reciprocal space using the cluster variation method was done by Sanchez in 
1982. EJ This k-space formulation of the CVM provides a significant improvement over the 
KCM formula, although a price is paid in increased computational effort. Thus, despite the 
improved treatment of fluctuations by the CVM, the method has been applied only to very 
few cases.li^ltJ An alternative approach, which we will refer to as the real space inverse CVM 
was introduced by Gratias and Cenedese.ll3 In this method, the effective pair interactions 
are obtained by fitting a few Warren-Cowley SRO parameters which, in turn, are related by 
a Fourier transform to the experimentally determined SRO intensity. A definite advantage 
of the real space inverse CVM is that the method involves only a modest computational 
effort, comparable to that of a regular CVM equilibrium calculation. On the other hand, 
the main drawback of the real space technique is that the fitting is done for only a few 
Warren-Cowley SRO parameters which, typically, are not enough to describe the measured 
intensities. To a great extent, this problem is solved by applying the inverse method in real 
space using Monte Carlo simulations, in which longer interaction ranges may be included 
without any fundamental difficulties.EJ'ta A problem of a slightly different nature is the fact 
that the determination of the Warren-Cowley SRO parameters from the experimental in- 
tensities, despite the fact that they are uniquely defined in terms of a Fourier transform, is 
not without errors. This step in the procedure, which of necessity requires the introduction 
of a real space cut-off in SRO, is an additional source of uncertainty on the derived effective 
pair interactions. 

Here we proposed a new method that carries out the fitting of the SRO intensity directly 
in k-space using the CVM formulation of the SRO intensity of Sanchez.^ Since the Warren- 



Cowley SRO parameters are not required, the k-space fitting circumvents two of tlie main 
problems of the real space method, namely the cut-off in SRO and the determination of 
the SRO parameters which, as mentioned, must be obtained by either fitting or by Fourier 
transformation of the experimental data. The method is applied to diffuse neutron scattering 
data obtained at several temperatures for a Fe — \Q.bat%Al single crystal^ A comparison 
between the real and k-space methods is also carried out for three different approximations 
of the CVM. 

The organization of the paper is as follows. In Section II we describe briefly the exper- 
imental aspects. The theoretical model and a review of the CVM formulation of the SRO 
intensity is given in Section III. Results are presented in Section IV with concluding remarks 
given in Section V. 

II. EXPERIMENTAL CONSIDERATIONS 

The experimental results, including standard corrections to the data, have been reported 
in previous publications.Q'Q Thus, here we briefly recall only the main experimental aspects, 
referring the reader to Refs. [ ^,|^ for further details. 

The diffuse scattering measurements were carried out at high temperatures on D7 at 
ILL (Grenoble, France), using a 10 mm in diameter and 20 mm long cylindrical single 
crystal of Fe — 19.bat%Al oriented along the [110] direction. Time of flight analysis was 
used in order to eliminate the phonon and critical magnetic contributions. Here we analyze 
four temperatures between 1273/ir and 15737^. The cross sections were calculated using 
the classical corrections (instrumental background, double scattering and absorption). The 
experimental points were averaged with a weight inversely proportional to the square of 
their error to obtain a regular grid of points (mesh of 0.1 reciprocal lattice unit step). 

In order to obtain the SRO intensity, the experimentally determined neutron scattering 
cross section was corrected for contributions due to static and dynamic atomic displace- 
ments. The correction for static displacements was carried out using the method based on 
the formulation of Borie and Sparks,B°l whereas dynamic attenuation was treated using the 
standard Debye- Waller theory. It should be noted that the accuracy of these corrections 
affect critically the reliability of the effective interactions derived from the diffuse intensity. 

In terms of the experimental scattering cross section, the corrected intensity in Laue 
units is written: 
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where N the total number of atoms and Fiaue is the usual normalization factor given by 
FLaue = c(l — c){bAi — bpeY , with c the aluminum concentration, and bAi and bpe the nuclear 
diffusion lengths of Al and Fe, respectively. In Eqn. (P, line is the incoherent scattering 
and B the Debye- Waller factor. 

Following Borie and Sparks,By the corrected intensity is approximated by a first order 
expansion in the atomic displacements Up, where the subindex p stands for the set of in- 
tegers {pi,P2,P3) and the location of lattice site p is given by Rp = -{pi,P2,P3)- In this 



2n 
approximation, the corrected intensity at the reciprocal space vector k = — (/ii, /12, h^) can 

a 
be written as: 

3 

Icorrik) = a{k) + Y^hiQiik) (2) 

— * — * 

where a{k) is the SRO contribution we are interested in, and the Qi{k) are related to the 
Fourier transform of the atomic displacements. We next consider each of these two terms 
separately. 

The SRO intensity a{k) is given by the Fourier transform of the Warren-Cowley SRO 
parameters a^. 

a{k) = X! '^p cos(A; ■ i?p) (3) 



with a-p defined by: 






where cTq and a-p are occupation numbers at the origin and at site p, respectively. These 
occupation numbers take values 1 or —1 if the lattice site is occupied by Fe or Al atoms, 
respectively, and the brackets <> stand for configurational averages. 

The quantities Q{k) = [Qi, Q2, Q3] in Eqn. (H) are given by the first order displacement 
parameters jp-. 

Q[k) = E 7p sin{k ■ Rp) (5) 

p 

In turn, the displacement parameters 7p are defined by: 



7p = -^ E E ^ P2,p(aa') < A«-' > (6) 

" <T=±1 a'=±l 



^ Laue 



where < A-Up'^ > is the average relative displacement between atoms of type a and a' 
(i.e. Fe or Al) separated by i?p, and where p2,p(ccr') is the probability of finding atoms 
a and a' separated by distance _Rp. We note that these probabilities are c^ + c(l — c)ap, 
(1 — c)^ + c(l — c)ap and c(l — c) — c(l — c)ap for AlAl, FeFe and FeAl pairs, respectively. 
The quantities a{k) and Q{k) are invariant to any symmetry operation of the space 
group of the reciprocal lattice, whereas the diffuse intensity Icorrik) is not. Consequently, 
the experimental points can be grouped into families of equivalent k-points having the same 
values of a{k) and Q{k), but with different intensities Icorr{k). Families having as many or 
more k-points than unknowns {a{k) and Q{k)) can be used to separate both contributions 

from Icorr{k). 

Since for the (110) reciprocal plane Qi = Q2, there are at most 3 unknowns on this plane. 
We note that at special points the number of unknowns is further reduced by symmetry to 
either 2, when Qi = Q2 = Q3, or 1, when Qi = Q2 = Q3 = 0. The maximum number 



of independent measurements possible at some of the equivalent points is 6. However, the 
measurement range in reciprocal space is limited by the neutron wavelength and by the lack 
of accuracy in the diffuse scattering around the Bragg peaks due to the mosaic structure of 
the single crystal and to different parasitic contributions. These limitations reduce strongly 
the occurrence of 6 equivalent measurement points. In fact, in the more frequent case there 
are 4 experimental points for three unknowns. In such cases (24 measured points), the 
unknowns were solved using a weighted least square fit, with a weight inversely proportional 
to the square of the experimental error. Thus, for each such point in the first Brillouin zone, 
one obtains a value a{k) with an error bar, Aa{k), and a measure of the quality of the fit, x^, 
given by the weighted sum of the square of the difference between the measured and fitted 
values. In order to take into account the quality of the fit, the error bar Aa{k) was multiplied 
by X when the latter was larger than one. For 18 k-points, the number of unknowns is equal 
to the number of experimental points. In these cases the system was solved and the error 
bar was multiplied by the maximum x of the first class of points. Finally, for ten points, 
mainly located around the origin and at the edges of the measured region, there are less 
measurements than unknowns and, therefore, it is not possible to separate a{k) and Q{k) 
from the measured intensity. 

A point of clarification concerning the Debye- Waller factors is also in order. These were 
measured from the attenuation of the fundamental Bragg peaks through neutron and gamma 
ray scattering.Ej We note, however, that the attenuation of the Bragg peaks is related to the 
sum of the dynamic, Bd(T), and static, Bs, Debye-Waller factors, i.e. Biro = BdiT) + Bg, 
whereas the attenuation of the short range order intensity is given by their difference, B = 
BdiT) — Bg. In order to separate these two contributions, the static part was assumed to 
be temperature independent. Since this approximation results in large errors in Bd{T) at 
high temperature, the restriction of a temperature independent static Debye-Waller factor 
was subsequently lifted by including a correction AB{T) in the SRO attenuation:Ea 

B = BdiT) -Bs + ABiT) (7) 

The correction ABiT) was determined by exploiting the symmetry of the corrected 
intensity as well as using a least square fit of the experimental data at equivalent k-points. 
The latter procedure was carried using points for which the number of unknowns is smaller 
than the number of measurements. 

The values of ABiT) obtained by the two different methods are very similar at high 
temperatures, indicating that the Borie and Sparks model is adequate to describe the scat- 
tering cross section. Departure are seen at low temperatures, likely due to the magnetic 
short range order contributions which are not explicitly considered in the analysis. Thus, 
the high temperature data reported in Refs.[ Q are considered to be the most accurate, and 
we will consequently restrict the present analysis to them. 

Figure 1 shows the experimental SRO intensity at 1273A' after all corrections have been 
made, together with the fit obtained with 8 Warren-Cowley SRO parameters (solid line). 
In what follows, we will refer to these quantities as the experimental Warren-Cowley SRO 
parameters, although they are actually the result of fitting the a(A;). In the fit, a value 
ao = 1.147 was obtained, which is to be compared with the exact value of 1. Each integer 
step on the horizontal axis of Fig. 1 corresponds to a scan along the (OO/is) direction, with 
hs between and 1, for fixed values of hi = h2- We note that the SRO data at the other 



temperatures analyzed here, namely 1373i^, 1473K and 1573ii', look qualitatively the same. 
Furthermore, at these temperatures, a similar quality of the fit is also achieved with 8 shells 
of a„'s. 

A notable feature of the SRO intensity is the small peak observed near the origin. It 
is not clear from the experimental data, however, whether this maximum is real or it is 
simply an artifact due to the cut-off imposed on the SRO parameters and/or to the fact 
that the fitting was carried out on only one plane in reciprocal space. The presence of 
this maximum, whether real or not, poses interesting questions on the use of the real space 
inversion methods. We will return to this point after discussing the theoretical models used 
to describe the SRO intensity in alloys. 

III. THEORY 

In this section we develop a simple statistical model in order to calculate the observed 
SRO diffuse intensity in terms of the effective pair interactions. We begin by describing the 
Hamiltonian, which is to be solved at finite temperatures using different levels of approxima- 
tion of the CVM, continue with a brief description of the approach used to calculate the SRO 
intensity, and conclude with a description of the inverse methods employed to determine the 
pair interactions. 

A. The Hamiltonian 

We describe the Fe—Al system with a simple Ising model in which the magnetic moments 
are localized on the Fe atoms: 



p,p 



+ \ E ^pV ^p^p' (8) 

p.p' 

where S^ = ±1 is the spin at site p. In Eqn. (|]), Vp^, and Jpp/ are, respectively, effective 
chemical and exchange interactions between sites p and p'. Averaging over the magnetic 
degrees of freedom, the interacting part of the alloy Hamiltonian can be written as: 

^2 = 2 E ^pp' ^p^p' (9) 

p,p' 

where the effective pair interaction is: 

Vpp' = ^pV - ^^pp'<5p5p'> (10) 

We note that in terms of interatomic potentials between different chemical species, the 
effective interactions are given by: 

Vpp' = \ [Vp^p^^^ + V^T - 2 V^T] (11) 



Thus, in the present model, the unknown magnetic SRO is included in the effective pair 
interactions, which now should depend on temperature.a Furthermore, since the product 
Jpp' < SpSpi > in Eqn. ( [I0| ) is alway positive, the Vpp/ are smaller than the chemical 
interactions V^p/- 

B. Cluster Variation Method 

The Hamiltonian in Eqn. (P) is solved at finite temperatures using the Cluster Variation 
Method. In this approximation, the alloy free energy is written in terms of the probabilities of 
finding different atomic configurations on clusters of lattice sites. A given cluster, consisting 
of n lattice points (pi, P2, •••, Pn), "will be characterized by a set of indices n, describing the 
geometry and the orientation of the cluster, and by its location in the lattice, fn. As needed, 
the three features of the cluster, type and orientation n and location fp, will be collectively 
denoted hj (3 = [n, rp]. Here we adopt the convention that the location of the cluster is 
given by its the center of gravity: 

r„ = ^ E Rp (12) 

"' pe/3 

For a given cluster (3 = [n, fn] there are 2^ atomic configurations, each one described by 
the set of occupation numbers a/3 = [o"p^, ■■■o-p^]. The probability of occurrence of any one 
of these configurations will be called p^{a^). 

With the energy given by Eqn. {^, the CVM free energy functional takes the form:cJ 

•^ = ^ J2 ^pp' < ^p^p' > + E ^p < ^p > 
^ P,p' p 

-keT E «/3 E Ppi^p) ln{pp{ap)) (13) 

where we have included a staggered effective chemical potential field pp, which will be used 
in the next subsection to describe general concentration fluctuations. The coefficients a^j in 
Eqn. ([131) are central to the CVM. After choosing a set of maximum clusters {7}, which 
determines the level of approximation, the coefficients a/3 are calculated from the following 
recursive relationsiBH 

E'«" = 1 (14) 

qD/3 

where the sum is restricted to clusters a contained in the set of maximum clusters 7, i.e. 
for clusters a C {7}. Furthermore, Eqn. (|l^ is valid for all clusters /3 included in the 
approximation, i.e. /? C {7}. The coefficients ap have the full symmetry of the disordered 
lattice and, therefore, they are independent of the location and orientation of the cluster. 
Furthermore, these coefficients are uniquely defined once the set of maximum clusters is 
chosen. 

The free energy functional, Eqn. (0), must be minimized with respect to the cluster 
probabilities Pp{(yp) subject to normalization constraints: 



E P«(^a) = l (15) 

(7a 

and, for all clusters /3 C a, subject to the self consistency requirement: 

Ppi^p) = E Pa(o'a) (16) 

where the sum is carried out over the occupation numbers dp = ±1 for all lattice sites p 
contained in a but outside cluster (3. 

A general description of the cluster probabilities, in which Eqns. (|I3]) and (jT^) are strictly 
obeyed, can be achieved by using a complete and orthogonal basis set of characteristic 
functions.cil'E3 These functions, in turn, can be used to describe any function of configuration, 
of which the probabilities P/si^is) are a particular case. For binary systems, the characteristic 
functions, which we will denote cr^, are defined for each cluster (3 in the lattice by the product 
of the occupation numbers CTp over all points p contained in /3:El 

^P=U ^P (17) 

pe/3 

In this basis, the cluster probabilities take a particularly simple form:cil 

1 



P/3l^/3J 



}(Ta<(^a> 



(18) 



where the multisite correlations, < (Tc >, given by the expectation value of the characteristic 
functions, form a set of linearly independent configurational variables. These multisite 
correlation can be conveniently used in the free energy minimization: 

'^ (19) 



d < as > 



P 

The minimum of JF defines the equilibrium free energy, F, as well as the equilibrium values 
of the multisite correlation functions, < cr^ >. The latter, of course, have the symmetry 
of the equilibrium phase. In the usual case in which the effective chemical potential field is 
spatially uniform, the number of correlation functions < <t^ > is relatively small. In such 
cases, the set of Eqns. ([T9| ) can be solved efficiently using fast numerical algorithms based, 
for example, in the Newton-Raphson method. In the next subsection, we will study the 
response of the system to small arbitrary variations in the effective chemical potential. The 
analysis leads to the fluctuation spectrum of the concentration variable < cXp >, which is 
directly related to the observed SRO intensity. 

C. Fluctuations 

We note that solution of Eqns. (|T^ yields a set of correlations or, equivalently, Warren- 
Cowley SRO parameters for only those pairs included in the set of maximum clusters {7}. 



Typically, the Warren-Cowley parameters obtained in this manner are not sufficient to re- 
produce the experimental SRO intensity. Thus, the required long-range pair correlations 
are calculated using a different approach, based on the ffuctuation spectrum obtained with 
the CVM free energy, Eqn. (|l^). For the sake of completeness, we summarize here the main 
results needed for the calculation of SRO intensities using the CVM, and refer the reader 
to the original work for further details.lij Following Ref . [ |14| , we write the pair correlations 
as: 

< o-pCTp/ > - < cTp >< o-p/ >= -ksT- — (20) 

5/ip9^p' 



A Fourier transform of Eqn. (^) yields:t3 



<\a{k)\^>-\<a{k)>\^ = -kBT ^ ^ ^ (21) 

(9/i(A;)9yu(— A;) 

with: 

<yik) = E ^P e-*^" (22) 



and 



^ E ^-'■''^4- (23) 



a/i(A;) V ^^p 



Actual calculation of the right hand side of Eqn. (|2T|) requires the Fourier transform of 



the matrix of second derivatives of the free energy with respect to the correlation functions: 

where the sum is over all the locations f^ and f^' in the lattice of clusters of type n and n' 
and where -Fn.n'l'^n — "^n') is given by: 

i^n,n'(r„ - f„0 = 7^ 7. (25) 

with 13 = [n, rj and /?' = [n', r^']- 

— > 

Since the effective chemical potential /i(/c) is the conjugate thermodynamic variable of 
the point correlation function < cr^k) >, the second derivatives in Eqn. (^) are given by:l3il 

NF,-Hk) (26) 



dii{k)dn{—k) 



where F-^ l{k) is the diagonal element corresponding to the point correlation function (n 
n' = 1) of the inverse of the matrix -Fn,n'(fc)- 

9 



Combining Eqns. i 
of the SRO intensity: 



an p6|), we obtain the desired result to be used in the fitting 



a{k) 



NkBTF.-ljk) 

I- < Go >2 



(27) 



The procedure leading to Eqn. (^) is equivalent to expanding the free energy around 
equilibrium up to second order in the correlation functions, and describing the fluctuation 
spectrum in the Gaussian approximation. As mentioned, Eqn. ( P7[ ) yields a significantly 
improved description of the SRO intensity over the Krivoglaz-Clapp-Moss formula. We note 
that although accuracy is generally gained as the cluster size increases, the size of the matrix 
of second derivatives Fn,n'(^) can rapidly grow beyond practical limits. In the next section, 
we will describe results for a CVM approximation that include clusters of up to nine points, 
which yields a complex second derivative matrix of size 1088 x 1088. 

D. Inverse k- Space Method 



The simplest procedure to obtain effective pair interactions from the measured SRO in- 
tensity consists in calculating the Warren-Cowley SRO parameters a-p from the measured 

— * _____ 

a{k), and then proceed to solve for the interactions using Eqns. (|TP[) with the proper con- 
straints on the pair probabilities. As mentioned, this real space inverse method has been 
applied to several cases using both the CVM and Monte Carlo simulations. Here we propose 
an alternative approach that does not require the computation of the real space Warren- 
Cowley SRO parameters but, instead, obtains the pair interactions by fitting a{k), calculated 
using Eqn. (pT]), to those obtained experimentally. The procedure consists in minimizing 
the following sum of errors squared: 



X 



1 

M 



M 

E 



amikj) - adki 
Aa{h) 



(28) 



with M the number of measurements, am{ki) and adki) the measured and calculated SRO 
intensities, respectively, and where Aa(/cj) is the experimental errors determined as described 
in Section II. 

In order to estimate the errors in the effective interactions Vn obtained from the least 
square fitting, we consider the probability distribution for the SRO intensity to be of the 
form: 



M 



P{{a{k)}) oc W exp 



i=l 



6a{ki 



Aa{ki 



(29) 



with 6a{ki) = a{ki)— < a{ki) >. 

The corresponding probability distribution for the effective interactions, up to second 
order in Vn, is then given by: 
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p({v;})ocn exJ-E^™,.^^ 

n [_ m,n ^Vm lWji 

with 5Vm = Vm— < Vm > and Am,n given by: 



(30) 



_^ 1 dajki) dajki) 

i=i Aa{ki) dVn dVm 

Thus, an estimate of the error in the effective interactions can be obtained from the eigen- 
values, Xn, and the matrix of eigenvectors, i'n,m,, of A^, „: 

AK = X E %^ (32) 



m V ^n 

where we have included the factor x, which is typically larger than one, in order to take into 
account the quality of the fit in the error of the effective interactions. 

IV. RESULTS 

The methods described in the previous section were applied to the three different approx- 
imations of the CVM, which will be referred to as the tetrahedron (T), the cube-octahedron 
(C-0) and the cube-rhombohedron-octahedron (C-R-0) approximations. The maximum 
cluster for each of these approximations are shown in Fig. 2. In the T approximation, only 
first and second neighbors are involved. The C-0 approximation, which has been used 
previously in the context of the real space inverse CVM,Q'Q extends the range up to fifth 
neighbors. In this approximation, however, the fourth neighbor pairs are excluded since 
they do not belong to the maximum clusters (the body-centered-cube and the octahedron). 
This difficulty is resolved by the C-R-0 approximation since the additional rhombohedron 
cluster includes fifth neighbor pairs. 

A. Real Space Method 

First, we apply the standard inverse CVM. For comparison purposes we have carried out 
the calculations in all three approximation, although the T approximation is not expected 
to be very reliable in this case since only two Warren-Cowley parameters are used in the 
fit. The results are shown in Table I where the first column lists the Warren-Cowley SRO 
parameters extracted from the experimental data for the first 5 neighbor shells. We recall 
that 8 shells of cui's have been used in the experimental analysis of the data.™ We also 
note that in the C-R-0 approximation all the a^'s listed in Table I are reproduced exactly, 
whereas in the C-0 approximation there is no information on a^ and in the T approximation 
only the first 2 a^'s are reproduced. 

The values of Vn are comparable in the three approximations, with the notable exception 
of V2 that is three to four times smaller in the C-R-0 than in the other two approximations. 
A direct test of the quality of the inverse CVM results can be obtained by calculating a{k), as 
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described in Section III.C, using the interactions of Table I. The results of these calculations 
for the three approximation are compared to the experimental data in Fig. 3. 

Overall, the quality of the fit for all three approximations is not very satisfactory. The T 
approximation (dotted line in Fig. 3), with only two effective interactions, gives a very poor 
agreement near the maximum at (0.5,0.5,0.5). On the other hand, the C-0 approximation 
(dashed line in Fig. 3) reproduces this maximum quite well, although it fails to reproduce 
the increase in intensity predicted by the experimental a„'s (see Fig. 1). As mentioned, the 
poor fit near the origin may not be particularly important since, based on the k-space data, 
the existence of a maximum near the origin is questionable. Nevertheless, this failure of the 
C-0 approximation to reproduce the diffuse intensity maximum illustrates two problematic 
aspects of the real space inverse CVM: i) The fact that the cut-off in SRO may introduce 
spurious effects in the SRO intensity and, ii) The fact that fitting just a few of the Warren- 
Cowley SRO parameters is not enough to reproduce all the features of the SRO intensity. 

The inverse CVM in the C-R-0 approximation (full line in Fig. 3) works very well around 
(0.5,0.5,0.5), reproduces the diffuse intensity maximum near the origin, but it gives poor 
agreement around (001). 

Finally, in order to illustrate the problems that may arise from using the C-0 approxi- 
mations, in which fourth neighbor pairs are not considered (i.e. V4 = 0), we have calculated 
the Warren-Cowley SRO parameters using the C-R-0 approximation with the interactions 
obtained in the C-0 approximation. The results are shown in the last column of Table I. 
We see that the 0:4 predicted by the interactions of the C-0 approximation differs in sign 
from the experimental one, given in the second column in Table I. 

B. Reciprocal Space Method 

The difficulties encountered in the application of the inverse CVM to the Fe — Al data 
are likely due to the uncertainty and errors in the experimental determination of the q;„'s. 
Although this difficulty can in principle be resolved by measuring more points on several 
reciprocal space planes, a more elegant solution is to fit the SRO intensity directly in k- 
space. The resulting pair interactions at 1273/^, together with the predicted a^'s, are 
shown in Table II. We see from the table that the interactions, particularly for the nearest- 
and next-nearest-neighbor pairs, are significantly different from those obtained using the 
real space method. Another noticeable difference is that the interactions shown in Table 
II change very little as we go from one approximation to the next. In Figure 4 we show a 
comparison between the results of the real space (open symbols) and reciprocal space (full 
symbols) methods for the T (triangle), C-0 (circles) and C-R-0 (squares) approximations. 
The lines connecting the V„'s obtained in the C-R-0 approximation are drawn as an aid to 
the eye only. 

The SRO intensity corresponding to the interactions of Table II are shown in Fig. 5. We 
see that the overall fit has improved considerably relative to that obtained with the real 
space method. Furthermore, the three approximations give essentially the same intensity. 
This is understood from the fact that the dominant interactions are the first and second 
neighbors (see Table II). Thus, the T approximation is already sufficient to capture the 
main features of the diffuse intensity. The increase in intensity near the origin predicted by 
the experimental q;„'s is not present in any of the three approximations. 
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Since the dominant effective interactions are between first and second neighbors, the 
SRO intensity for the other temperatures investigated here were analyzed using the T ap- 
proximation. The results for Vi and V2 as a function of temperature are shown in Fig. 6. 
There is a gradual increase in the effective interactions with temperature which is expected 
from the decrease in magnetic SRO (see Eqn. ([T0|)). 

V. CONCLUSIONS 

The determination of effective interactions from measured SRO intensities requires con- 
siderable manipulation of the experimental data. In this process, the guidance of reliable 
theoretical models is invaluable. In particular, we have shown that the ability to calculate 
accurately the diffuse intensity in reciprocal space allows the determination of the effective 
interactions without the intermediate step of obtaining the Warren-Cowley SRO parameters 
in real space. Consequently, the method proposed here avoids the need to introduce a real 
space cut-off in SRO. 

The study of fluctuations in reciprocal space by means of the CVM allows the descrip- 
tion of statistical correlations beyond the range of effective interactions. This establishes a 
fundamental difference with the KCM formula, derived from the Bragg- Williams approxi- 
mation, in which the fluctuation spectrum is determined entirely by the Fourier transform of 
the interaction potentials. Furthermore, these statistical correlations, which extend beyond 
the range of the effective interactions, are describe well using relatively small clusters. For 
example, for the Fe — Al studied here, and due to the fact that first and second pair inter- 
actions are dominant, we find that the Tetrahedron approximation suffices for an accurate 
description of the experimental data. 

As an illustration of the advantages of the k-space versus the real space methods, we 
note that the k-space method in three different approximations produces essentially the same 
results in terms of the effective potentials and predicted SRO intensities. In contrast, the 
results of the real space CVM are, for the same data, significantly different to each other and 
to the results of the k-space analysis. Particularly revealing is the fact that increasing the 
number of Warren-Cowley SRO parameters used as input in the real space inverse CVM does 
not necessarily produce better results. For example the inaccuracies of the T approximation 
which, incidentally, works well in the k-space analysis, can in principle be explained by 
the fact that only ai and 02 are considered. Increasing the q;„'s to five in the C-R-0 
approximation improves the results in some regions in reciprocal space but also introduce 
large errors in other regions and, possibly, spurious effects such as the increase in intensity 
near the origin. In fact, for the case studied here, the C-0 approximation, that neglects 
the fourth neighbor interactions altogether, appears to produce the more sensible results 
in terms of the predicted intensities. However, a closer look shows that the interactions 
obtained in the C-0 approximation give the wrong sign of the fourth neighbor Warren- 
Cowley parameter. 

Finally, we point out that the k-space CVM analysis of the SRO intensity calls for a 
different experimental approach. Whereas the real space methods require measurements 
of a{k) in many points in the irreducible Brillouin zone in order to perform the Fourier 
transformation, only a few points are actually needed when using the k-space CVM analysis. 
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Thus, a more efficient strategy to obtain effective interactions is to measure the diffuse 
scattering at fewer points in reciprocal space but with higher accuracy. 
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FIG. 1. Experimental short-range order diffuse intensity (symbols) and the fit (solid line) ob- 
tained using 8 shells of Warren-Cowley short-range parameters a„. 
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FIG. 2. Clusters used in the three different approximations of the Cluster Variation Method: 
Tetrahedron (T) approximation (points 1-2-5-7); Cube-Octahedron (C-0) approximation com- 
bining the body-centered cube (points 1-2-3-4-5-9-10-11-12) and octahedron (points 2-3-4-5-6-9) 
clusters; and the Cube-Rhombohedron-Octahedron (C-R-0) approximation combining, in addition 
to the two clusters in the Cube-Octahedron approximation, the rhombohedron labeled 1 through 
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FIG. 3. Short range order diffuse intensity obtained using the real space inverse cluster variation 
method in three approximations: Tetrahedron (dotted line); Cube-Octahedron (dashed line) and 
Cube-Rhombohedron-Octahedron (full line). The symbols are experimental points. 



FIG. 4. Effective pair interactions as function of distance, in units of the lattice constant a, 
obtained using the real space (open symbols) and the k-space (full symbols) methods described 
in the text. For each method, the cluster variation approximations used were the Tetrahedron 
(triangles), Cube-Octahedron (circles) and Cube-Rhombohedron-Octahedron (squares). 
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FIG. 5. Short range order diffuse intensity obtained using the reciprocal space cluster variation 
method described in the text in three approximations: Tetrahedron (dotted line); Cube-Octahedron 
(dashed line) and Cube-Rhombohedron-Octahedron (full line). The symbols are experimental 
points. 
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FIG. 6. Temperature variation of the effective pair interactions for first (Vi) and second (V2) 
neighbors obtained using the reciprocal space method in the Tetrahedron approximation. 
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TABLES 

TABLE L Pair interactions Vn for Fel9.5at%Al at 1273i^ obtained using real space inversion 
for three approximations of the cluster variation method. Errors in the estimated pair interactions 
are shown in parenthesis. Experimental Warren-Cowley short range order parameters are given in 
the second column. Those in the last column were calculated in the C-R-0 approximation with 
the interactions of the C-0 approximation. 





Exp. 


T 


C-O 


C-R-O 


n 


Ui 


Vn{meV) 


Vn{meV) 


Vn{meV) 


ai 


1 


-0.0943 


19.7(0.2) 


18.5(0.3) 


18.5(0.2) 


-0.0944 


2 


0.0114 


4.2(0.2) 


3.6(0.2) 


1.1(0.2) 


0.0113 


3 


0.0345 




-2.4(0.1) 


-4.8(0.1) 


0.0351 


4 


0.0097 






-3.5(0.1) 


-0.0091 


5 


-0.0020 




2.0(0.1) 


2.0(0.1) 


-0.0019 



TABLE II. Pair interactions Vn and the corresponding Warren-Cowley short-range order pa- 
rameters for Fel9.5at%^/ at 12TiK obtained by fitting the SRO intensity in k-space using three 
approximations of the cluster variation method. Errors in the estimated pair interactions are shown 
in parenthesis. 





T 




C- 


-O 




C-R- 


-0 




n 


Vn{meV) 


an 


Vn{meV) 




"n 


Vn{meV) 




an 


1 


28.7(1.9) 


-0.1158 


25.7(4.7) 




-0.1107 


23.7(1.9) 




-0.1067 


2 


10.8(1.7) 


-0.0044 


8.7(3.7) 




-0.0005 


7.1(0.4) 




0.0047 


3 






-0.7(1.2) 




0.0326 


-0.9(0.1) 




0.0332 


4 












0.3(0.1) 




-0.0080 


5 






0.6(1.8) 




0.0076 


1.6(0.1) 




0.0016 
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